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Abstract 

To calculate the temperature at which a first-order cosmological phase tran- 
sition occurs, one must calculate F^{T), the free energy of a critical bubble config- 
uration. F^(T) is often approximated by the classical energy plus an integral over 
the bubble of the effective potential; one must choose a method for calculating the 
effective potential when l^" < 0. We test different effective potential approxima- 
tions at one loop. The agreement is best if one pulls a factor of fi^/T^ into the 
decay rate prefactor [where /z^ = and takes the real part of the effective 

potential in the region V^" < 0. We perform a similar analysis on the 1-dimensional 
kink. 
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1. Introduction 



Thermal Tunneling and the Critical Bubble Free Energy 

A scalar field theory whose potential V has two local minima may tunnel out of the 
false vacuum (0^) by the nucleation and subsequent growth of bubbles of true vacuum 
While we will refer to V as the "classical" potential, it may arise in part from integrating out 
other particles in the theory, e.g. gauge bosons'^', so V may have implicit temperature (T) 
dependence. The nucleation rate per unit volume in the static limit {RT ^ 1) is calculated 
in the Gaussian approximation (i.e. to 1-loop order) to be'^"^' 



V ~ V ^ 2 



-1/2 



e-^=/' (11) 
det[a2 + /x2] V ■ ; 

where //^ = is the classical energy of the critical bubble, a static and spheri- 

cally symmetric field configuration 0(r), of radius R, which extremizes the classical action'^' 
subject to periodic boundary conditions in Euclidean time. The determinants range over 
a complete basis of fluctuations about the classical solution (0(r) or subject to the 
same boundary conditions, < is the eigenvalue of the "breathing" mode about 0(r). 
The second term on the RHS of Eq. (1.1) is from Affleckl^l, and the | is from analytically 
continuing the breathing mode integration!^! . 
With the periodic boundary conditions, 

aet[d' + V"m = exp I f; In [{27rnTf + a;|] 1 , (1.2) 

tn=-oo j ) 

where the cj| are eigenvalues of [— + V"{(j))], and the (cUj*)^ are eigenvalues of [— + /x^]. 
We use the identity!^' 

I ^ In [{2nnTf + a;'] = f + rin(l - e"'^/^) + C = Tin [2 sinh (^)] + C . (1.3) 

n 

The constants C cancel out in Eq. (1.1). The uj_ contribution is then traditionally pulled 
back into the prefactor. The 3 "translation" modes {n = and Uq = 0) are not treated 
correctly above; they actually give V{E^/27rT)^^'^ in the prefactor^^], and the remaining lOq 
contribution (from n^O modes) gives in the prefactor. This gives 

r ^ /^y/^ \u;_\/2T .../^ 
V 2tt\2ttTJ sin(|u;_|/2T) 

where the "traditional" bubble free energy 

ptrad = £;^ + AF*;"^ = + AFf "'^ + AF^^"'^ , (1.5) 
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= E'^V^ + ^^t""' = E'^ln "ZD ■ (1-6) 

J ^ Vl-e / 

Primes on the sums in Eq. (1.6) indicate omission of the translation and breathing modes 
{ujj, j = 1 — 4). Counterterms F^* are discussed below. 
We now define^ 

psub = E^ + AFll\, = F^ + AFi""^ + AF|"\ (1.7) 

^psub ^ ^p^trad^ ^p^uh ^ ^pU;ad _ 4Tln{T/n) . (1.8) 

Now Eq. (1.4) becomes 

V 27r V27rry sin(|u;_|/2r) ' ^ ' ' 

We will find that the effective potential approximation most closely approximates F^"''. 

The Effective Potential 

The sums in Eq. (1.6) are often approximated by treating the fluctuations locally as 
plane waves to get an effective potential V^^rp = Vi+Vrp, then integrating [V^i+t(0)~^i+t(0/)] 
over all space. No attempt is made to remove the 4 translation and breathing modes. In 
Eq. (1.6) one substitutes 

E-/d^x^^^, uj^^^k^ + V'i^), cuj^^/k^, (1.10) 
and one finds, with rn^ = V"{(f)), 

rri4 roo 

VA<I>) = ^ I{m/T), Hy)^J^ dxxMn(l-e-V-^+^^j . (1.12) 
The expansion of I{y) for real y < 2% is^^'^l 

,2 \ f^' 



4 2 4 

45 12^ 6^ 32 



,f 4(2fe)!C(2fe + l) /VY 



;i.i3) 



I^This is somewhat like removing the lowest 4 cj^'s from the sums in Eq. (1.6), in addition to the lowest 4 
uj/s, since their contribution to F*™"^ is -4[^ + Tln(l - e"**/^)] w 4Tln(T//i). 
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where C3 = | + 21n(47r) — 27 5.4076. We choose a renormahzation scheme in which all 
divergent graphs are precisely cancelled by counterterms so that at zero external momenta, 
— ^1 (0/) — — (^^d there is no wavefunction renormahzation), 1^1 specifically: 



net 



-1 

64^ 



d^x { [4A^ + i/x"] + m2[4A2 - 2//^] + 



2 -In 



4A2 



m^=V" 



(1.14) 



In the region < 0, we must modify these results to give a real answer. For Vi we will 
always take the real part of Eq. (1.11). For let us keep the first equation of Eq. (1.12), 
but replace I{m/T) by /("^s)(|m|/T) where 



/(-^)(y) ^^-^Y' + Y' [a + Mn(y^)] - ^ [in(y^) -c, + c\ + 



(1.15) 



Methods we consider are then parametrized by {a,b,c}. The most common and obvious 
method (A) is to take the real part of Eq. (1.13), corresponding to {a = 6 = c = 0}. 
Another method (B), proposed in ref. [9], replaces the lower limit of integration in Eq. (1.10) 
by A; = Im{m} (eliminating fluctuations with wavelengths longer than the bubble thickness). 



and corresponds to {a = | — | ln(2), b 



6' 



0}. 



The Derivative Expansion 

For configurations 0(x) which vary slowly, the effective potential approximation is the 
leading term in a derivative expansion of the free energy. The next term (at high T) is'^^'^^^ 



T 



1927r 



d^xm"W^I 



m 



(1.16) 



and again we take the real part (Method A) when necessary. More terms are given explicitly 
in ref. [11]; they become increasingly divergent at = 0, where the derivation breaks down 
(because an integration by parts becomes invalid). Also, no attempt is made to omit modes. 
The usefulness of Eq. (1.16) is thus highly suspect, but we note that derivative corrections 
are predicted to be C(T^). 

Scales, Approximations, and Goals 

Our generic tree-level potential will be quartic in with 0^ = 0, V"{0) — jj?, and 
(t)^ — a. By rescahngt^^l = (70, x — x/ji, and T — jiT, we can rewrite the 4-action Sq as 



Sr. 



T 



1 

T 



d^x 




1^2 _ ^^^3 ^^^4 

2^ 3^ 2^ 



(1.17) 
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K > 1 is a dimensionless parameter; >1 (degenerate minima) is the thin-wall limit, while 
larger k gives thicker bubbles. With tildes indicating dimensionless results, 

= (aV/x) E^, AF,^^ = fiAF,^^ . (1.18) 

The loop expansion^^'^l is an expansion in (/x/cr)^ and T. It is sometimes claimed that higher 
loops should eliminate the complex terms in F^°^, but this cannot be generally true since 
the higher-loop contributions are suppressed by these arbitrary parameters. Henceforth we 
will drop the tildes and work in the rescaled theory {i.e. set /i = o" = 1). 

We always use the static approximation'^^' {RT ^ 1) and the 1-loop approximation. 
In Section 3 we will use the thin- wall approximation, R ^ 1. At times we will make high- 
temperature expansions, requiring T > 1 (note the thin-wall and high-temperature limits 
together imply the static limit). We are examining the validity of the effective potential 
approximation. 

In this paper we will study several systems: the 1-dimensional kink, the thin-wall bubble, 
and two thick-wall bubbles. We will calculate AF^ and AF^ for each system exactly [F^^^ in 
Eq. (1.8)], in the effective potential approximation [F*'°* from Eqs. (1.11-1.12), using different 
methods to calculate /(^^f) in Eq. (1.15)], and using the next term of the derivative expansion 
^pder (1.16)]. 



2. The 1-Dimensional Kink 

Classical Results 

We warm up by calculating the free energy of a kink in 1 spatial dimension: I^^l 

g = vi<p) = i<p\i-<pr. (2.1) 

The potential is that of Eq. (1.17) with k, — 1. The kink solution is (up to an arbitrary shift 
in coordinate) 

^{x) = |[1 - tanh(ia;)], V"{^{x)) = 1 - fsech^^a;) . (2.2) 
Eq. (2.1) allows us to convert integrals over x into integrals over 0: 

For example, the classical energy is 

Note that in ID [compare to Eq. (1.18)] E^ — jJt.a'^ E^ and AF^ j, — /lAF-^ j., so with scales 
restored El^ — iJ,a^/6. 
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Method 




Tln(r) 


T 


ln(T) 


1 






sub 


-.2138 





-1.2425 


.4775 


1.0522 





.0036 


pot(A) 


-.0916 





-2.1145 


.4775 


1.0522 





.0036 


der(A) 


-.0916 





-2.1730 


.4775 


1.0522 





.0036 


pot(B) 


-.0916 


0.4495 


-1.7222 


.4775 


1.0522 


.0045 


.0036 



Table 1: Kink free energy in low- and high-T regimes. 

Exact Results from the Eigenvalue Sum 

The solutions to the eigenvalue equations (setting /jl — 1) are known:'^^'^^' 



= ^ , k,^ ^^'''^ , S{k) = 27r - 2 tan-i {k) - 2 tan'^ {2k) , (2.5) 

where we have imposed vanishing boundary conditions on a box of length L, so s is a positive 
integer. We drop the translation mode eigenvalue cui, there is no negative eigenvalue in ID. 
In the continuum limit, 

^X + X Vd^^— "2^ + ^ ' 
^p^rad ^T\n(l- + J^" ^^Tln (l - e"^^/^) . (2.6) 

In our renormalization scheme the ID counterterms analogous to Eq. (1.14) are 
F^* = ^ y dx { [4A2 + 1] + [2 + 2 ln(4A2)] - m^} [^Z^ = ^ P + 6 ln(4A2)] . (2.7) 

(This differs from ref. [15] by S/Sn due to different renormalization schemes; also note their 
= /x72.) We define AF^^"^ = AFf '^'^ and AF|"^ = AF*r«'^ - Tln(T//x), and find 

^psub ^ J_ _ A ^ -.2138 . (2.8) 
' 4^3 Stt ^ ^ 

AF-^ = -(In ^/l2) T + A ln(r) + + T"^ + • • • (2.9) 

where = 1 + 21n(47r) - 2-f ^ 4.9076, and C(3) ^ 1.2021. These results are in the row 
marked "sub" of Table 1. 
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ID Effective Potential and Derivative Expansion Results 



The ID effective potential for real m is'^' 

Vi-^Hm') + ^^, V^^-I{m/T) (2.10) 

For < we replace I{m/T) by I^'^^^\\m\/T) where 

j(neg)^y^^ ^^j,^^^y2^) _ f^yS-j _ ^ ^ . . . ^2.12) 

6 L J 8 

Method A gives {a = 6 = c = 0}, and Method B gives {d = 1 - ln(2), 6 = "i, c = 0}. 

We integrate (the real part of) Vi from Eq. (2.10) over all space, using Eq. (2.3), to 
get AFf ^^^-^ = —.0916, which differs significantly from AFf"* = —.2138 (note each result is 
renormalization-depcndcnt, but the difference is not). This difference, which was calculated 
in ref. [15], dominates the low-T regime. 

A similar integral for the high-T expansion gives 

AFf;f ) = ln[2(V3- V2)^] T + ^ ln(T) + + ^ T'^ + • • • , (2.13) 

ZTT 87r OZTT"^ 

as shown in the line marked "pot (A)" of Table 1. Note that the difference between the true 
result and the potential approximation no longer lies in the constant term, but only (as far 
as we have taken the expansion) in the T term! It is 

AFlfj, - AFf^*^^^ = - In 4V3(v^- v^)^] T = .8720 T . (2.14) 

The next term of the derivative expansion [analogous to Eq. (1.16)] is 

AFf " - AFI"* = ^ /" dxm-^ V'(m') = ^ ln(v^ - V2)T = -.0585 7 (2.15) 
96 J 48 

as incorporated in the third line of Table 1. It is a very poor approximation to Eq. (2.14)! 

Results from Method B are given in the fourth line of Table 1; these are also unsatisfac- 
tory. In fact, the choice {a — 1.940, 6 = c = 0} in Eq. (2.12) would give the correct ("sub") 
results, but it is not clear if there is any physics in this choice. 
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3. The Thin- Wall Critical Bubble 



Classical Results 



For K close to (but larger than) unity in Eq. (1.17), the solution to 

= V'{4>) (3.1) 
is a thin-wall bubble, given approximately by the kink solution in the radial coordinate, 
Eq. (2.2) with x — r — R and R » l.'^l The tree-level critical bubble energy has volume and 
surface terms: 



An J r^dr ^ ^ 



-^7TR^\V{l)\+i7TR^E^ 



2 iT-lD 



(3.2) 



where E^^ ~ | was given in Eq. (2.4), and |T^(1)| = (k — l)/6. We extremize to find the 
bubble radius R and energy E^, 



Stt 



27ri?2 



The wall thickness is 0{1) (i.e. lJi~^). It can also be shown^^] that u;^ ^ — 2/i?^, so the static 
and thin- wall limits imply that the third factor of Eq. (1.4) is near unity. 

Exact Results for a Domain Wall 

In the thin-wall limit, the surface frcc-cncrgy density fi^ = AF^ j,/ (47ri?^) of the bubble 
wall equals that of a planar domain wall^^"'!. We can thus solve the eigenvalue equation in 
Cartesian coordinates, using Eq. (2.5) for the radial wavenumber k^, and plane waves for the 
tangential A;^, to get 

2 2 27r 



n 



JO 

+ 
-1 



2ti 



5(v/A^) 



-4 



fsub 

Jt 



327r2 Vv^ 
1^ dk^ 
27r 



TT \k^ + l 4A;2 + 1 
6 ) = -.02474 , 



^0 



In [1 - 6-^=*/^] + In 



1 _ g-V^/^ 



+ 



^0 



87r2 327r2 



dk^ 

TT 



ln(4A2) 



X 



+ 



k^ + 1 Ak^ + 1 



In 



1 _ g-Vfef+fcr + l/T' 



(3.4) 

We have performed the /j, integral numerically, and fit to an expansion in T~^; the results 
are shown in Table 2 in the row marked "sub"."'" 



■i- These results are also useful for the study of second-order phase transitions, in which the domain wall free 
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f 
I\ 








Jl+T 








Method 






Tln(T) 


T 


ln(T) 


1 






sub 


-.02474 


-1/4 





.15215 


-.01900 


-.03712 





-.00012 


pot (A) 


-.00661 


-1/4 





.15452 


-.01900 


-.05612 





-.00012 


der(A) 


-.00661 


-1/4 





.15187 


-.01900 


-.05612 





-.00012 


pot(B) 


-.00661 


-1/4 


.00864 


.16409 


-.01900 


-.05612 


.00006 


-.00012 



Table 2: Thin-wall bubble free energy density for low- and high-T. 

Effective Potential and Derivative Expansion Results 

Results from integrating the effective potential, and the next term of the derivative 
expansion, over the bubble [again using Eq. (2.3)] are shown in the rest of Table 2. Using 
the general of Eq. (1.15) gives 

/f°*y ^-\T^ - (.0518 6)Tln(T) + (.1545 + .0259 a - .0242 6) T 

- (.0190) ln(r) + (-.05612 - .000514 c) ' ^^'^^ 

Matching this to the true /^"^ gives the coefficients {a, b, c} shown in the first line {k = 1) 
of Table 3.§ 

We see "derivative corrections" are 0{T). The derivative expansion prediction, ff^j, 
from Eq. (1.16), is a reasonable approximation to them in this case. 



K 


a 


b 


c 


1 


-.0913 





-36.974 


1.5 


.2834 





- 1.424 


2.5 


.4188 





- 0.180 



Table 3: /("^») parameters that make AFf^j. = AF^fj 



l+T' 



energy density is set to zeioM"^^ Restoring units 

^2 rp2 



f[<t> 



wall] 



Y~~f + -^^^^^ l^'^c - -01900 HTjfJ') 







•. That is, the critical temperature is a bit higher than the 



giving, for <C tr, = ^2/3a + 0.3/x 
leading result which is in the literature. 

^First subtracting the derivative correction of Eq. (1.16) from Ai^f"^ would give a values of .0109, .3877, 
and .5128, respectively. For the kink it gives a = 2.070. These results are no more enlightening. 
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A F 








A F 






Method 




T 


2 


Tln(T) 


T 


ln(T) 


1 


sub 


-2.13 


-78 


.61 





49.52 


-5.193 


-15.64 


pot(A) 


-2.65 


-78 


.61 





45.47 


-5.193 


-16.12 


der(A) 


-2.65 


-78 


.61 





43.98 


-5.193 


-16.12 


pot(B) 


-2.65 


-78 


.61 


4.76 


49.73 


-5.193 


-16.12 


thin-1 


-1.81 


-78 


.61 





47.84 


-5.974 


-17.65 


thin-2 


-4.97 


-50 


.27 





30.59 


-3.820 


- 7.46 



Table 4: Thick- wall bubble free energy for k — 1.5. 
4. Thick- Wall Critical Bubbles 

Classical Results 

From Eq. (1.17), the (scaled) potential (Fig. 1) is 



3^2^ 



(4.1) 



Larger k > 1 gives thicker bubbles. The minima are at = and = 1, with V"{0) = 1 
and V"{1) = 2k — 1. The bubble profile is the solution to 



(f)" + 2(f)' /r = 0(1 - 0)(1 - 2K(f)) . 



(4.2) 



Fig. 2 plots 0(r) and V"{r) for k = 1.5 and k = 2.5. From ref. [12], the classical energy is 
approximately 



4.85a 



a ( 2.4 .26 

^4 V ^1^^ (l-a)2 



9fi; 



a 



[1 + 2kY ■ 



(4.3) 



Exact, Effective Potential, and Derivative Expansion Results 

Our method of calculating the exact free energy F^*"*, formally given by Eq. (1.7), is 
described in ref. [8]. The results for k — 1.5 are in Table 4, and for k — 2.5 in Table 5**, along 
with effective potential and derivative expansion approximations. Thin-wall predictions are 
also shown for two values of R: one chosen to give the correct coefficient ("thin-1"), and 
one given by Eq. (3.3) ("thin-2''). Finally, the parameters in /("^f) needed to match the 
effective potential approximation to the exact result are given in Table 3. 



**In our fit to the data, we allowed a T ^ term, not shown, and constrained the T^, Tln(T), and ln(T) 
terms. 
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Method 




^2 


Tln(r) 


T 


ln(T) 


1 


sub 


-1.34 


-24.90 





17.17 


-1.408 


-4.60 


pot(A) 


-1.009 


-24.90 





14.05 


-1.408 


-4.64 


der(A) 


-1.009 


-24.90 





13.35 


-1.408 


-4.64 


pot(B) 


-1.009 


-24.90 


2.48 


15.60 


-1.408 


-4.64 


thin-1 


-0.572 


-24.90 





15.15 


-1.892 


-5.59 


thin-2 


-0.553 


- 5.59 





3.40 


-0.424 


-0.83 



Table 5: Thick- wall bubble free energy for K — 2.h. 

5. Conclusions: A New Prefactor, and Derivative Corrections 

We have tested the effective potential approximation to the critical bubble free energy. 
The agreement is best if one pulls a factor of /x^/T^ into the decay rate prefactor, Eq. (1.9), 
and takes the real part of the effective potential in the region V" < (Method A). That is, 
ppot{A) f^Yoselj approximates F^"^^ = F^^"''^ — 4T'ln(T//i). Table 3 shows that no single set 
of /("^f) parameters {a, 6, c} does consistently better than Method A. With scales restored, 
= 0{a'^/n), AF^fj, = 0{T'^/n), and "derivative corrections" are 

AFf^^ - AFf;^^) = 0{T) . (5.1) 

This difference is numerically fairly small, and very poorly predicted by the derivative ex- 
pansion [Eq. (1.16)]. In summary. 



T _ a' f Y'^ |^J/2r 



X^JB,] J:-.",':^. e^^-"VT , (5.2) 



V 2n\2nTj sin(|u;_|/2T) 

where X is a dimcnsionlcss number representing derivative corrections, typically 10~^ to 10^. 

In ID, where AFf"^ is only 0{T), derivative corrections [still 0{T), and numerically 
larger] are much more significant than in 3D. 
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Figure 1: The potential V(0) for several /c's. 



Figure 2: Thick-wall bubble profiles 0(r) and V"{r). 
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